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Abstract 

We investigate the general geometric properties of the surface of 
infinite red-shift corresponding to the event horizon of the static and 
axisymmetric solution of the Einstein vacuum equations that only pos- 
sesses mass M and quadrupole moment Q. The deformation of the 
Schwarzschild surface r = 2M produced by the quadrupole moment 
is shown, and the range of values of this multipole moment is speci- 
fied, which preserves a regular, closed, continuous and differentiable 
surface. Some thermodynamic consequences and speculations ensuing 
from our results are discussed. 
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1 Introduction 



As it is well known [1], the only static and asymptotically-flat vacuum space- 
time possessing a regular horizon is the Schwarzschild solution. All the other 
Weyl exterior solutions [2]-[6], exhibit singularities in the physical compo- 
nents of the Riemann tensor at the horizon. 

For not particularly intense gravitational fields and small fluctuations, 
deviations from spherical symmetry may be described as perturbations of 
the Schwarzschild exact solution [7]. 

However, such perturbative scheme will eventually fail in regions close 
to the horizon (although strictly speaking the term "horizon" refers to the 
spherically symmetric case, we shall use it when considering the r = 2M 
surface, in the case of small deviations from sphericity). Indeed, as we ap- 
proach the horizon, any finite perturbation of the Schwarzschild space-time 
becomes fundamentally different from the corresponding exact solution rep- 
resenting the quasi-spherical space-time, even if the latter is characterized 
by parameters whose values are arbitrarily close to those corresponding to 
Schwarzschild metric [8]-[T2]. This, of course, is just an expression of the 
Israel theorem (for observational differences between black holes and naked 
singularities see [13], and references therein). 

Therefore, for strong gravitational fields, no matter how small the mul- 
tipole moments of the source are (those higher than monopole), there exists 
a bifurcation between the perturbed Schwarzschild metric and all the other 
Weyl metrics (in the case of gravitational perturbations). 

Examples of such a bifurcation have been brought out for the 7 metric 
and for the space-time possessing only monopole and quadrupole 
moments (hereafter M-Q space-time [23|, [25], [2S]), in the study of the 
trajectories of test particles for orbits close to r = 2M [27], [28] . 

The influence of the quadrupole moment on the motion of test particles 
within the context of Erez-Rosen metric [29] has been investigated by many 
authors (see [30]-[33] and references therein). 

The purpose of this paper is to study further the properties of the M-Q 
space-time in order to bring out the bifurcation mentioned above, and to 
contrast our results with those obtained for other Weyl metrics. 

The rationale for the choice of the M-Q space-time is based on the fact 
that its relativistic multipole structure (particularly that of a sub-class of this 
solution M-Q*^^^ [21], which, being an exact solution, corresponds to the first 
order in quadrupole moment in M-Q) may be interpreted as a quadrupole 
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correction to the Schwarzschild space-time, and therefore represents a good 
candidate among known Weyl solutions, to describe small deviations from 
spherical symmetry. 

In this paper we shall see that unlike in the 7-metric [23], the area does 
not vanish at the horizon. However, the surface gravity diverges as is the case 
for the 7-metric [31]. We shall calculate explicitly the surface of infinite red- 
shift corresponding of the M-Q space-time and determine the range of values 
of the quadrupole moment for which such a surface is close, continuous and 
different iable. Next we shall calculate the surface gravity which exhibits a 
singularity along the symmetry axis as we approach the event horizon. Some 
thermodynamic consequences derived from these results are discussed. 



In [35] a so-called Multipole Symmetry Adapted coordinate system (MSA) 
was defined which allows us to write the metric component of pure multi- 
pole solutions in General Relativity, in a form resembling the classical multi- 
pole potential. The existence of such a type of coordinate system was proved 
in [35] . The relation between the MSA system of coordinates (r, y = cos 9) 
and the Weyl standard coordinates [R, u = cos 0) is given by the following 
expression: 



fn{^), fi'n(i^) being polynomials in the angular variable u explicitly calculated 
in [35J, where the consistency of this expansion, and the convergence of the 
series is discussed. In particular, for the case of the solution with a finite 
number of multipole moments given by the monopole and the quadrupole 
moments, the corresponding MSA coordinates can be written as follows: 



2 The M-Q solution 



r = R l + ^/„(a;) 



71=1 




(1) 



n=l 



r 



= r, [1 - qP2iys)\l - ^i^^yt + ^yl - 9)At- 

+ (^(105y,^ - myl + 1) - q\l2yt - Qy] + 1)) + . . . 
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l-(l-ys)(?A^ + |(7y.^ + 9)At+ 



+ ^A^ ((175y2 + 57) + g(126 - 336y2)) + . . . 



(2) 



where Pn{y) stands for the Legendre polynomials, r^, and ys denote the 
standard Schwarzschild coordinates, As = M/vs, and q = Q/M^, {Q stands 
for the quadrupole moment). 

The relevance and the physical and mathematical interest of these sys- 
tems of coordinates become evident from the fact that they are related with 
the existence of certain symmetries of the static and axisymmetric vacuum 
field equations, which generalize the action of those symmetry groups on the 
classical equations [36] • For the monopole solution (spherical symmetry) the 
MSA system of coordinates is exactly the standard Schwarzschild coordi- 
nate, and the associated symmetry group of the field equations allow us to 
determine univocally the system of coordinates by means of solving certain 
Cauchy problem. In a recent work |37|| the MSA system of coordinates has 
been used to calculate relativistic corrections of the perihelion advance for 
orbital test particles. 

The M-Q solution of the static and axisymmetric vacuum equations (see 
P?] . |24] ) can be written in the corresponding MSA system of coordinates 
X = (r, y) as follows 



9rr (•^) 



■1 + 2 
1 



1 - 2A 

1 M2 

M2 



1 + {l-2\) J2^^'U,{y,Q) 

i=3 

oo 

i + ^VA(y,Q) 

oo 

l + ^VT,(y,Q) 



A 



i=3 



(3) 



where A = M/r, and f/j, Di and Tj denote polynomials in the angular variable 
y of even order depending on the quadrupole moment Q. The explicit form 
of some polynomials are given in 
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3 The event horizon 



In this work we want to make use of the MSA system of coordinates to 
calculate the event horizon of the exact monopole-quadrupole solution of the 
static and axisjTumetric vacuum Einstein equations. 

In MSA coordinates the gu metric component of the M-Q solution takes 
the form 

goo = ~l + 2Vir,y), (4) 

with V{r,y) = h ^P2(y), which is a reminiscence of the classical mul- 

tipole potential but now constructed with the relativistic monopole mo- 
ment M and the quadrupole moment Q respectively. The event horizon 
of this static solution is given by the equipotential surface V{r, y) = (3 with 
P = 1/2. Equivalently, the surface of infinite red-shift defined by the condi- 
tion Qu = ^°'C,a = {^°' being the time-like Killing vector) can be calculated 
by solving the following algebraic equation 

-2Mr'^ -2QP2{y) = 0. (5) 

Therefore, the appropriate root of this equation is a function r = r{y) which 
represents the generatrix curve of revolution surface describing the event 
horizon of the axisymmetric solution. Let us note that this curve is symmetric 
with respect to the equatorial plane due to the symmetry of the solution and 
therefore, the study of r = r{y) can be restricted to the value of y E (0, 1), 
or equivalently the angular coordinate 9 G (0,7r/2). 

In order to determine the real roots of this third degree algebraic equation 
we must define the different domains of the equation in terms of the values 
for y = cos 9. It depends on the value of the quantity 

d = qP2{y) [qP2{y) + + 2«)' (6) 

with the notation a = qP2{y), a = 8/27, and q = Q/M^, as follows: 

DI : d<0,ae{-2a,0), 
DII : > 0,a G (-oo,-2a) U (0,oo), 
Dili : d = 0,a = -2a,a = 0. (7) 

The region defined by DI leads to three real roots of the algebraic equa- 
tion ([5]), as well as the region Dili (where at least two of these roots are 
equals), whereas region DII provides only one real root. 
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Let us first analyze the region Dili: the solutions of the algebraic equa- 
tion for this region are the following 



rrf=o(y) = ri(y) = 2M 
rd=o{y) = r2,3{y) = 2M 



+ (a + a 



.1/3 



,1/3 



(8) 



It is easy to show that we have positive real roots leading to a good matching 
with the region DII only if a = 0, which corresponds to the angular value 
= 1/ a/3 : the other possible point of this region is 1/2 



which corresponds to the condition a ■ 
However, the roots for that value are: 



2a, as can be seen from Figure 1. 



ri(Z/2°) 



2M 



^^2,3 (2/2 



AM 



(9) 



and they only match with the unique real root for the region DII [d > 0) 
providing a negative value of the curve: 



rd>oiy) = M 



{a + a + Vdy^^ + {a + a- Vd)^ 



/3 



•rd>o{y 



2M 



(10) 

Therefore we must limit the value of the parameter q within the range 
—2a < q < Aa, where the value of 1/2 is not real, to preserve the continuity 
of the curve r = r{y). 

Next, as indicated in ([7]) , the region DII {d > 0) splits into two different 
open intervals given by the conditions a > and a < —2a. The condition 
a > leads to y < y\ (for g < 0) or ?/ > (for g > 0) whereas the 
condition a < —2a cannot be fulfilled if the regularity range of the parameter 
q e {—2a, 4a] holds. 

Hence, we can conclude that the different domains for the evaluation of 
the real roots of (jSD are given by 



q G (0,4a] 



q E (-2a, 0) 




(11) 



(12) 



y>yi 
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Figure 1: 




— m 

The function d = d{a) in terms of a 
Legendre polynomial of degree 2 in ® • The two points where d vanishes 
the variable y = cosO. are a = and a = -2a {a = 8/27). 

Finally, the region DI {d < 0) has the following three real roots: 
r{y) = -M(^2cos- + lj = rd<o, 

r{y) = ^M(2cos^±y3sin^ + l) =r±, (13) 

where k is defined by cosk = 1 H — . Let us note that these solutions r-t, 

a 

TfK^Q are well defined only if —2a < a < which in fact agrees with the range 
of a for this region ([7]). Additionally k = arccos ^1 H — j exhibits a good 

behavior since it is uniquely valued at the boundary point y = y^, i.e. k has 
a single value (k = 0) at the limit a — )■ [y = yl) since k G [0, tt] where 
-2a < a < 0. 

The curve r = r{y) matches appropriately at the boundary between the 
regions DII and Dili where 

rd>o(y = y?) = 2Af = ri(y = y?) (14) 

and the curve is completely continuous for all range < ?/ < 1 if we consider 
the curve r = riy) to be equal to r^^o ffT3|) at the region DI (rf < 0), the 
other roots r± being neglected. 



8 



In conclusion, M-Q solution possesses an event horizon defined by a gen- 
eratrix curve which is continuous for the whole range < < tt as follows: 



gG(0,4a]: r{y) 



ge(-2a,0): r{y) 



rd<o 
rd>o 

rd>o 
rd<o 



0<y<y'i 
yl<y<i 

< y < 
yl<y<i 



(15) 



(16) 



Let us note that the continuous curve r = r{y) must be a piecewise-defined 
function since k is not defined for the region d > DII where a > 0. 
Nevertheless it can be shown that functions rd<o and rrf>o are equivalent 
since we can rewrite r^j^o dlS]) by taking into account that 



a K 
cos K = 1 H — ^ cos — 
a 3 



2ai/3 L 



+ a + ^Y'^ + {a + a-^Y'^] (17) 



and therefore 

rd<a{y) 



M 



+ a + y/df'^ + (a + a - ^Y'^ + 



(18) 



which represents the root rd>o- 

Finally, we can prove that the surface of the event horizon defined by 
the generatrix curve (IT5l)-( lT6|) is class C^, i.e. r(|/) is different iable and the 
derivative is continuous. The derivative of the generatrix curve (fT5|) -( IT6|) with 
respect to the angular parameter 9 is given by the following expressions: 



dr, 



d<0 



dO 
drd>o 



—AM qy . k 
y^ : sin — 

^ 3 3 



1 - y' 



;qMy 



+ a + VdY'^ - (a + a - Vdf/"^] , (19) 



de ' ' Vd 

and the continuity of the derivative at the boundary {y = y^) is fulfilled: 



dr, 



d<0 



dy 
dy 



y=yl 



y=y'i 



4M qy k 
lim : sm — 

'^^o 3 3 



3gM 



lim 



qMy 



2^3 

+ a + Vdy/'' ■ 



[a + a 



Vdy 



/3 



3qM 
2^3' 

(20) 
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In addition, for any value of the parameter q within the range considered 
(— 2«,4a], the derivative of the curve r = r{y) is also continuous at the 
points intersecting the symmetry axis y = ±1, as well as at the equatorial 
plane {y = 0), since both derivatives of ra^o and rd<o are zero at those points. 

Hence, the surface of the event horizon is continuous and differentiable ev- 
erywhere for the range of values of the quadrupole moment Q G M^{—2a, 4a]. 
Figures 2 to 4, show different aspects of the generatrix curve as well as the 
surface of the event horizon. 

It should be observed that although pathological behaviour might be 
expected for large (absolute) values of the quadrupole moment, it is not 
clear why such a behaviour appears outside the range established above. We 
currently do not know if there is any physical significance behind that range. 



The area of the event-horizon surface r = r{y) for the M-Q solution is given 
by the following expression: 



and grr, gyy, gifiip are the metric components of the M-Q solution given by 
expressions ([3]) in the corresponding MSA system of coordinates. As we have 
previously seen, the event horizon for the non-spherical case depends on the 
angular variable 6 (cos 6' = y) for the range defined in f|T5|) -f fT6|) and hence its 
derivative with respect to the variable y is different from zero. Nevertheless, 
if we expand the expressions of the surface ffT^ - ffTB]) as well as its derivative 
( !T9|) in terms of the quadrupole parameter q, we obtain 



4 The Area 




(21) 



where 




(22) 




(23) 
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and we could approximate the event horizon to be equal to the Schwarzschild 

surface r ~ r,, = 2M and — ~ 0, since we are interested in the behaviour of 

dy 

the area for very slight deviations from spherical symmetry. Equivalently we 

4 K 

can say that r{y) = 2M + e, where e = -M(— 1 + cos — ) representing the dif- 

3 3 

ference between the event horizon of the M-Q solution and the corresponding 
surface in the spherical CcLSG, IS cL small quantity of order q'^. 

Therefore the area of the event horizon could be defined as follows: 

Ac^2n J y/gyyg^^dy, (24) 

Since we want to know how different is the area of the event horizon for 
the non-spherical case for slight deviations with respect to the Schwarzschild 
space-time, let us consider the M-Q^^^ solution. 

M-Q(^) is a subclass of the M-Q solution, that represents a small de- 
formation of Schwarzschild for a small parameter q. In fact, the M-Q^^^ 
solution corresponds to the first order in the parameter q in an expansion of 
the solution M-Q in power series of that parameter. 

When written in prolate spheroidal coordinates {x, y}, the event horizon 
is defined by x = 1 (r = = 2M) if q < 8/5, and metric functions of the 
solution are: 



where 



A = l + -g -1 + 



(3?/2- l)(3a;2 



4" V 4 
5 / 4_ 

— y^ 
15 2, 



C = -—qx{l-y 



D = -'4^{i-y')-l9(^' + y')j^,^, (26) 



T 

^It can be seen from figures 3-4 that — = 2 + p where p = t/M with —2/3 < p < 1/3, 
and p ~ ior q — leading to r{y) — 2M. 
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and the terms of order higher than q have been neglected. 
The area of the event horizon is 



A=\mvM^ / e-^"^ Ve27(x2 - y^){x^ - l)dyd^, (27) 

Jq J-l 



or, using (125])- (|26 



A = 27rM^ hm /' {x + iy-^''/^~^\x - l)i+(^/2-A)gB+D/2^^_ ^28) 



Then, taking into account that 



hm(a; + l)i-(^/2-A) ^4 ^ hm(x - = 1, (29) 



we have that 



A ~ 2n4M^ j e~l'^^^+^y''>dy, (30) 



which evaluated up to order q produces 



A~47r(4M')(l-^g) + 0(g2). (31) 



In the spherically symmetric case ( 1311) brings us back to the well known 
result 

Asch. = levrM^. (32) 

5 The surface gravity 

The surface gravity k may be defined from the time-translation Killing vector 
C as 

= -I (V,e.) (V-D • (33) 

In a static, asymptotically fiat space-time it is the acceleration of a static 
observer near the horizon, as measured by a static observer at infinity. 

Also it may be defined more generally for any surface outside the horizon 
in terms of the four-velocity of a fiducial observer at rest with respect to the 
frame for which our metric functions are defined |38]. At event horizon of 
course both definitions coincide. 
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For the case of a static and axisymmetric line element (diagonal matrix) 
the surface gravity k is given by the following expression: 



(34) 

where the index i = 1,2 sums for the non-axial coordinates. 
By taking into account the expression (jl]) we have that: 

de goo = —^6 cos 6 sin 9, (35) 

and therefore the surface gravity for the M-Q solution at the event horizon 
is given by the following expression 



MJl-2X-2X^qP2{y) 



(l-2A)(l + 3gA3p2(y))^ , 9\Wa-y') ^'^ 



l + (l-2A)f/ 1 + ^ 

(36) 

where A = M/r, and f/, D denote the series appearing at the metric functions 
grr and gyy respectively which are given by expressions (E]). 

In the spherical case we have g = 0, D = and f/ = producing 

which at the event horizon r = 2M, (A = 1/2) leads to the well-known 
expression for the surface gravity of Schwarzschild space-time: 

1 

^Sch = -TT7- (38) 

Let us now calculate the surface gravity for the M-Q solution by con- 
sidering the event horizon r{y) as an small deviation with respect to the 

r 

Schwarzschild surface r = 2M as follows: — = 2 + » where p = e/M, and 
A ~ -(1 - |) + O(p^). Hence, the expression (!36|) just at order zero in power 
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series of the parameter p leads to a divergent term aX y = = 1/ -\/3 since 
the surface gravity is given by 



An analogous conclusion follows for the M-Q^^^ solution. Indeed, in this 
case we have from fl25p - (|2^ 

5-00 = -e , 

— 1 

= M^e-*e^-^^, (40) 
from which it follows at once, 

where the following notation has been introduced 
X = A,., In - 5. 



X + 1 J — 1 

= Aln(|^)-B„. (42) 

In the spherical case = 0), as we approach the event horizon (x — )■ 1) we 
have: 

2 

x'^ — 1 

producing with ( l34l) 



X = -^ , y = 0, (43) 



1 



ks^^ = . (44) 

If we now consider the parameter q to be arbitrarily small but non- 
vanishing, then from fHTj) we obtain 

..2 4 15 [31/2 - 1 x(3t/2 - 1) /x - 1 

~ \ q h -^^ ^In' 

(x^ — 1)^ — 1 



x2-l ■ 2 ""Vx+l 
+ 0(g2 



1 /_ 4 

^ 0{q^), (45) 
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and therefore the surface gravity for the M-Q*^^) solution is given by the 
following expression 



where 



4M 



A = 



1 + -g(l - 15?/2) + A 



-q lim(x 

2 X— >i 



(46) 



(47) 



The important point is that everywhere outside the symmetry axis [y ^ ±1), 
A vanishes (at the event horizon), whereas for the points on the symmetry 
axis 1/ = ±1 at the horizon, A goes to infinity thereby producing a divergent 
surface gravity. This particular behaviour exhibited in the limit of (H71) is 
characteristic of the M-Q space-time, and has been brought out before in 
[28]. 

It should be stressed that this divergence of the surface gravity on the 
event horizon at y = ±1, A is not a "directional" singularity, i.e. surface 
gravity would diverge at the points ?/ = ±1, on the horizon, no matter how 
we approach it, as can be seen from P7|) . Therefore such a divergence is not 
related to a breakdown of the coordinate system as seems to be the case for 
"truly" directional singularities [39] . 

Finally, it should be observed that consistency with third law of thermo- 
dynamics would require surface gravity to go to zero when the area goes to 
zero. This condition is not satisfied for the gamma metric (see discussion in 
[33]). In our case however the area of the horizon remains finite. 

We shall next discuss about the consequences emerging from the results 
presented so far. 



6 Conclusions 

We have presented a detailed description of the event horizon for the M~Q 
solution. The range of values of q for which that surface is regular, closed, 
continuous and differentiable has been clearly established. 

We have next calculated the area of the event horizon. This quantity plays 
a central role in the definition of different bounds for black hole entropy (see 
[38] . [lQ]-[ll] and references therein). 

Now, if the gravitational field of the collapsing body, including small de- 
viations from spherical symmetry, is described by means of an exact solution 
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to Einstein field equations (instead of perturbations of Schwarzschild space- 
time) then it appears that in some particular cases (e.g. the 7 space-time) the 
area surface of the event horizon vanishes leading to the conclusion that as 
the body falls through the horizon, information about the collapsing body is 
stripped away, thereby resolving the information loss paradox |23]. However 
as we have seen before this is not the case for the M-Q solution, therefore 
the conclusion in [23] seems to lacks of universality, and does not apply to 
any deviation from the spherically symmetric case. Thus all the discussion 
about entropy bounds for spherically symmetric black holes remains basically 
unchanged for slightly distorted space-times (small q) in the case of the M-Q 
solution. 

On the other hand however, the result obtained for the surface gravity 
indicating that it diverges, points in the same direction that the one obtained 
for the 7 metric ^34j (nevertheless observe that in our case k only diverges 
on the axis of symmetry whereas in the 7 metric case it diverges everywhere 
except on the symmetry axis). 

Now, the consequences of such divergence might be important and de- 
serves a deeper analysis. 

On the one hand, the divergence of k would imply that the "inertial" term 
in the transport equation would grow unlimited as the horizon is approached. 

Indeed, as shown by Tolman many years ago [15] according to special 
relativity all forms of energy have inertia, and of course this should also 
apply to heat. Therefore, because of the equivalence principle, there should 
be also some weight associated to heat, and one should expect that thermal 
energy tends to displace to regions of lower gravitational potential. This in 
turn implies that the condition of thermal equilibrium in the presence of a 
gravitational field must change with respect to its form in absence of gravity. 
Thus a temperature gradient is necessary in thermal equilibrium in order 
to prevent the flow of heat from regions of higher to lower gravitational 
potential. This result was confirmed some years later by Eckart [16] and 
Landau and Lifshitz [17] . In the transport equation derived by these authors, 
the "inertial" term deduced by Tolman appears explicitly and is essentially 
the surface gravity multiplied by the temperature. It should be observed that 
exactly the same term appears also in transport equations derived from, more 
physically reasonable, causal thermodynamic theories (see [H] for details). In 
the case of Eckart and Landau and Lifshitz approach, the transport equation 
reads 

= -K/i"^(T;3 + Ta^), (48) 
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where q'^ is the heat flow four-vector, is the four-velocity of a fiducial 
observer mentioned above, a'^ is the four-acceleration associated to those 
fiducial observers, h^^" is the projector onto the three space orthogonal to 
V^, K denotes the thermal conductivity, and T denotes temperature. 

Therefore the unlimited growing of the "inertial" term as the surface of 
the object is closer to the event horizon would lead to an increasing inwardly 
directed heat flow vector, which in turn would produce a substantial increase 
of temperature. Now, it has been shown (see [H], [19] and references therein) 
that, as the system leaves the equilibrium, the inertial mass decreases by a 
factor (1 — a) where a is defined by 

kT 

a = , (49) 

where r is the thermal relaxation time and [i and denote the energy 
density and the radial pressure of the fiuid respectively. 

Obviously for this effect to be of some relevance in the dynamics of the 
system, evolution should proceeds in such a way that a approaches the critical 
value of 1. 

Now in c.g.s. units (omitting the pressure term which is always smaller 
than the energy density) we have 

r/i 81 [r] [^] ^ ' 

where [n] , [T] , [r] and [/i] denote the numerical value of these quantities in 
ergs~^ cm~^ K~^, K, s and gcm~^ respectively. 

Thus in order for a to attain values close to the unity we need extremely 
high values of thermal conductivity and/or temperature. In the past [5D] it 
has been suggested that neutrino radiative heat conduction might produce 
values of [k] as large as 10^^, which together with the expected values of 
temperature [T] ^ 10^^ at the last stages of massive star evolution, would 
lead to a ~ 1. 

Here we see that perturbations of spherical symmetry described by means 
of the M-Q space-time (this also applies to the 7 metric case), could provide 
a mechanism for attaining large values of temperature due to the increasing 
in the heat fiow as the source approach the horizon. 

On the other hand, starting from the generalized version of the surface 
gravity, it appears that its divergence as we approach the horizon, implies 
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the divergence of the Tolman mass (see [5H])- This resuh whose dynamical 
imphcations cannot be overemphasize, is also obtained for the 7 metric (see 
[5T] for a detailed discussion on this point). 

Before concluding, the following comment is in order: 

According to the black hole has no-hair theorem, in the process of con- 
traction all (radiatable) multipole moments are radiated away [52J . Therefore 
it could be asked what the interest might be to study singular horizons as 
the one corresponding to the M-Q solution? 

Nevertheless, the situation is more complex than it looks at first sight. 
Indeed, let us admit that in the process of collapse, all (radiatable) multipole 
moments are radiated away. Obviously, the mechanism of radiation, as any 
physical process, must act at some time scale (say Tmech.)- Now, if Tmech. is 
smaller than the time scale for any physical process occurring on the object 
(say Tphys.), then the appearance of a regular horizon proceeds safely. 

However, let us suppose for a moment that there is a physical process 
whose Tphys^ is of the order of magnitude of (or still worse, smaller than) 
Tmech.- In this case any physical experiment based on such process "will see" 
a singularity as the boundary of the object crosses the horizon, due to the 
always present fluctuations. 

Indeed, the fact remains that perturbations of spherical symmetry take 
place all along the evolution of the object. Thus, even if it is true that close 
to the horizon, any of these perturbations is radiated away, it is likewise 
true that this is a continuous process. Then, as soon as a "hair" is radiated 
away, a new perturbation appears which will be later radiated and so on. 
Therefore, since "hairs" are radiated away at some finite time scale, then 
at that time scale {Tmech.) there will be always a fluctuation acting on the 
system. 

Thus, unless one can prove that indeed Tmech. is smaller than Tphys. for 
any physical process, one should take into account the possible consequences 
derived from the presence of fluctuations of spherical symmetry (close to the 
horizon), in this later case such deviations from spherical symmetry might 
be described (at least in some cases) by the M-Q space-time. 
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Figure 2: Generatrix curve r = r{y), except for a factor 1/M, in polar 
coordinates, for different positive values of q. The vertical axis Z represents 
the symmetry axis y = ±1, whereas the horizontal abscissa axis X represents 
anyone direction on the equatorial plane orthogonal to Z. The curve for each 
value of q is drawn in different color depending on the corresponding domain 
up or down the boundary value y = yl, where all different curves are crossing 
themselves. 
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Figure 3: The whole curve r{y)/M for all values of the angular polar coor- 
dinate and different values of g (g = 0, g = 1.08518 and q = —0.53518). The 
curve with points represents the event horizon for the spherical symmetry 
case. 




Figure 4: The surface of the event horizon for different values of q. 
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